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The Lamb shift in muonic hydrogen was measured some time ago to a high accuracy. The 
theoretical prediction of this value is very sensitive to the proton-finite-size effects. The proton 
radius extracted from muonic hydrogen is in contradiction with the results extracted from elastic 
electron-proton scattering. That creates a certain problem for the interpretation of the results from 
the muonic hydrogen Lamb shift. For the latter we need also to take into account the two-photon- 
exchange contribution with the proton finite size involved. The only way to describe it relies on the 
data from the scattering, which may produce an internal inconsistency of theory. 

Recently the leading proton-finite-size contribution to the two-photon exchange was found within 
the external field approximation. The recoil part of the two-photon-exchange has not been con¬ 
sidered. Here we revisit calculation of the external-field part and take the recoil correction to the 
finite-size effects into account. 
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I. INTRODUCTION 

A few years ago the Lamb shift (2pr/2 ~ 2si/2) in 
muonic hydrogen was experimentally determined to a 
great accuracy [I|. The related theoretical prediction 
is much less accurate, suffering from a bad knowledge 
of the proton charge radius. Therefore, the experimen¬ 
tal data and theoretical expression have been utilized to 
determine an accurate value of the proton radius. Ac¬ 
tually, theory (see, e.g., [H-Q) disagrees with experiment 
and therefore the muonic-hydrogen value of the proton 
charge radius is inconsistent with results of hydrogen 
spectroscopy Q and evaluations of the elastic electron- 
proton scattering (see, e.g., il). 

The theoretical prediction for the Lamb shift in muonic 
hydrogen is basically of the form 

AAi = Cl -f C2 , (1) 

where the first term, dominated by the Uehling correc¬ 
tion and other QED effects, does not depend on the pro¬ 
ton rms electric charge radius Re, and the second term 
is dominated by the leading proton-finite-size correction. 
Meanwhile, there is a higher-order proton-finite-size cor¬ 
rection, which does depend on the distribution of the 
proton charge (or rather we prefer to discuss the electric 
charge form factor of the proton G_e( 9^)), but does not 
have the required form of const x R^. 
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The general problem in its calculation is that to find 
it one has to use a certain description of the charge dis¬ 
tribution (or GEiq^))- There is no satisfactory ab initio 
model. So, one either could apply a certain variety of 
models, which are not really related to experimental data 
and arbitrary to a certain extent, and to estimate some¬ 
how the uncertainty, or, alternatively, one could take ad¬ 
vantage of the existing experimental data. The first op¬ 
tion does not really work for the proton. As for the sec¬ 
ond option, the experimental data by themselves are not 
accurate enough and one cannot use the data directly. 
Instead one has to apply a fit. Unfortunately, and that 
is a real problem with the proton radius, the fits also 
produce their value of and those values from empiric 
fits do not agree with the value from the muonic hydro¬ 
gen (see, e.g., parameters of low momentum expansion of 
some empiric fits in Table HI. 

We have to clarify that the correction we discuss here 
is a small one, while the contradiction is huge. The very 
fact and the scale of the discrepancy does not depend on 
the way we treat the higher-order proton-finite-size term. 
Meanwhile, the eventual value of the muonic Re does 
depend on the way, we treat the higher-order finite-size 
correction, at the level of the uncertainty of Re- Here 
we revisit the correction. 

The higher-order proton-finite-size correction can be 
understood as a certain moment of the charge distribu¬ 
tion. Since the very mean radius of that distribution is 
under question because the results from the scattering 
and the muonic hydrogen are different, we cannot di¬ 
rectly apply the scattering results for these corrections 
(see, e.g., [Ml)- 
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Thus, the next-to-leading proton-finite size corrections 
are sensitive to detail of the charge distribution. A self- 
consistent consideration of those shape-sensitive contri¬ 
butions was performed in [l^ . The consideration was 
done within the external field approximation, which de¬ 
livers the leading part of the shape-sensitive contribu¬ 
tions. Here, we present a complete evaluation, including 
recoil effects. Those involve not only the charge distribu¬ 
tion, but also a distribution of the magnetic moment. 

There are two main ingredients of our consideration. 
We have to explain the method of a treatment of the 
higher-order finite-size correction within a theory of 
muonic hydrogen, which would be consistent both with 
the experimental data on the electron-proton scattering 
and with the muonic value of the proton radius. The 
method has been applied previously [l^ within the ex¬ 
ternal field approximation. Here we are to expand it to 
accommodate the recoil finite-size corrections, which we 
also need to consider. For that we follow [Ij, [T|. In 
both cases we start following cited papers, but modify 
the approaches. 


formalism, such as the effective Dirac equation (EDE) 
[3, [H, [13 • Within such an approach the TPE diagrams 
appear with certain subtractions (related to proton pole 
contributions) and on top of those subtracted two-body 
TPE diagrams one has to add a perturbative series of an 
effective one-body theory, e.g., the EDE theory with a 
proton on mass shell. The latter is presented by the Friar 
contribution, while the former is related to two types of 
the effects. Those include elastic (with subtractions) and 
inelastic (polarizability) parts (see Figs. [5] and [3]), which 
at a certain stage are to be considered together. 



FIG. 2: Two-photon-exchange: elastic contribution. Neces¬ 
sary subtractions, which follow from an EDE consideration, 
are assumed. The circles are for the vertex with the appro¬ 
priate form factors. 


II. GENERAL EXPRESSION FOR THE 
ELASTIC TWO-PHOTON-EXCHANGE 
CONTRIBUTION 


Let’s start with the expression for the proton-finite- 
size contribution including the recoil effects. They are 
described by the two-photon exchange (TPE). The lead¬ 
ing part of the higher-order proton-finite-size term is the 
so-called Eriar term. That is a result derived within the 
external field approximation (see Eig. [T|). The result for 
the correction was first derived in the coordinate ^ace 
[IM3 and next in the momentum space [Tsj l (see [3 for 
details of the derivation) 


AEpr{nl) = 
Ipr = 


16{Za)^ 


j Sio 

■ 


dq 


3 ’ 


( 2 ) 


(GEiq^)) -I-2G'E{0)i^ 


Here and throughout the paper we apply the relativistic 
units, in which h = c = 1, m stands for the muon mass, 
M is the proton mass and nir is the reduced mass. The 
momentum q is the Euclidean one. 
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FIG. 1: The diagram for the Friar correction to Lamb shift in 
muonic hydrogen. The appropriate subtractions are assumed. 


If we like to consider the TPE correction with all re¬ 
coil effects included, we have to use a certain two-body 



FIG. 3: Two-photon-exchange: inelastic contribution (the 
proton polarizability). 

Indeed, once we have well-defined elastic and inelas¬ 
tic quantities, we can work on their calculation indepen¬ 
dently. However, their definition is not a straightforward 
issue. As an example one can consider a calculation of 
the elastic part. It requires an appropriate on-shell ver¬ 
tex for the proton. One may deal with the Dirac and 
Pauli form factors Fi _2 or with the Sachs ones Ge,m- 
The difference between those two parametrizations leads 
to a contribution, which vanishes for the on-shell vertex, 
but is not equal to zero for the TPE diagrams in Eig. [21 
which involve the vertexes off-shell. In other words, the 
elastic part is not well-defined by itself. 

Meanwhile, both elastic and inelastic contributions 
have their specifics and their calculation is very differ¬ 
ent. Indeed, there are distinguished contributions into 
each of them and those contributions can be studied in¬ 
dependently of the definitions. One of them is the Friar 
term, the dominant part of the elastic TPE and we in¬ 
deed can study it independently of anything else. (As 
we mention, the ambiguity in a possible definition of the 
elastic term concerns usually only recoil effects, while the 
Eriar term is the result within the external field approx¬ 
imation.) Indeed, the actual expression for the elastic 
part for the two-photon-exchange (eTPE) correction is 
more complicated than the Friar contribution ([2]). Once 
we intend to complete the calculation of the elastic part 
and to add recoil corrections to the Friar term, we have 
to check the definitions. 
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Therefore, we have to explicitly quote the definitions 
of the elastic and inelastic TPE and to be sure they are 
consistent. While calculations of the elastic and inelastic 
parts of TPE have been done in a single paper for a few 
occasions [2ll - l^ . the first complete evaluation of TPE, 
where consistency of their definitions was discussed, ap¬ 
peared some time later in [l3j |. Unfortunately, it did not 
manage to produce consistent definitions of the elastic 
and inelastic part. The consistency was achieved later in 

M- 

The elastic part of the TPE contribution consists of 
a QED result for a point-like proton and of a proton- 
structure correction. We acknowledge that the point¬ 
like part has been already calculated and it is subtracted 
while calculating the eTPE term. The QED TPE term is 
the well-known Salpeter contribution (see, e.g.. 

It is for the proton without any anomalous magnetic mo¬ 
ment, which satisfies the conditions Ei(0) = 1,^2 (0) = 0 
or Gij(O) = Gm(0) = 1. 

The analytic expression for the eTPE contribution was 
considered in [3 [11 nn, m. Here, we re-arrange the 
expression for our purpose. The general structure is of 
the form [Il,[3 (cf- Ellll) 


AEeTPE(«0 = ~ 


16{Za)^mf Sio 

- JeTPE^ ) 


.^eTPE — .^Fr T Irec ; 

“ dq 


Ipr = 
Uec — 
.^E = 

Im = 
If = 


(GE(g")) -1-2G^(0)(72 

0 9 

.^E + + .fp , 


/o 9 

/o 


fM{m,M;q'^) (Gm(<7^))^ - 1 


^ ^fF{m,M;q'^) {Fi{q^)f - 1 (3) 


The Dirac form factor 


Fi{q^) = 


GE{q^) 


4mH 


Gm(<?') 


1 + Sd 


is more fundamental than the Sachs electric and magnetic 
form factors Ge and Gm- However, just the latter are 
used to be fitted in empiric-fit papers. 

The last term. Ip, was excluded in [3] from the elastic 
part and attributed to the polarizability. Later on, it was 
re-established by [3 as a part of the elastic TPE. Most 
of the existing results on the polarizability (3; HH) 1343 
are consistent with this definition^ or, as in case of [3| . 
can be re-adjusted to it. 


^ Actually, the very calculation of the polarizability in [T3l is con¬ 
sistent with the definition of the elastic part including this term 
as was shown in H- 


The first term apparently is the Eriar term, which is 
a non-recoil contribution. The other terms are for the 
recoil. As it happens in the two-body systems (see, e.g., 
various corrections to the Lamb shift in i) the external- 
field contributions and the recoil one are treated differ¬ 
ently and they have different subtractions. The structure 
G^ — 1 is a result of the subtraction of the expression for 
the point-like proton. 

The explicit presentation of the / factors [3 is 


fEim,M;q^) = 


fM{m,M\q^) = 


fF{m,M-q'^) = 


q 


(l2{Tp) l2{Tp)\ 1 


2(M — m) \ 

1 1 
/tZ ,/r, 


1 + Fb 


qTp 


2(M-m) 




ll{Tp) 7l(T-p) 


M / J 

72 (Tp) 72(rp)\ 1 


Tp J 1 + Tp 


M + m 


Fpll{Tp) 


( 4 ) 


where 


4M2 

„2 


71M 
72 ( 1 -) 


4m2 ’ 

(1 — 2r) (Vl + T — i/t) -I- a/t 

(l+T)3/2-r3/2-^y/7. 


( 5 ) 


We are going to treat electric and magnetic form fac¬ 
tors somewhat different and we are to decompose the 
Dirac form factor 


- 1 


(1+^^)2(G| - 1) + J^^^iGMGE - 1) 

+ (1 + p^)2(G'L - 1) (6) 


into terms of the Sachs form factors. We also note that 
the eTPE contribution takes into account two structure 
effects, which are a non-zero value of the anomalous mag¬ 
netic moment of the proton ^ 2 ( 0 ) = Gm(0) — 1 = k ~ 
1.792 85 and the finite size of the proton. It is helpful to 
disentangle them because the former can be found with a 
very hig h accuracy, while the latter as we will show (see 
also [l2|) needs for its calculation additional experimen¬ 
tal data and is not free from uncertainties. 
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Eventually, we arrive at 


Irp.r. — 


/« 

-^EF 

-fMl 

Im2 

where 


Ik + ^EF + ^M1 + Im2 

rOO rJ 

^ / “{(2 + ^)/mi + /M2} ) 

Jo 9 

^ -\fEF{m,M\q^) {Gsiq^)) -1 


/o 9 
dq 


/ —/mi { GM { q ^)) — (1 + 

9o 9 L 

f -j/m2 [Gm(9^)Ge(9^) ~ (1 + «^)] ) (7) 

9o 9 


/ef = /e + 


(1 + ^p)^ 


/e 


/mi = /m + 


/m2 = 


2to 


(l + 'TP)^ 


(1 + 'Tp)^ 
/e • 


/e 


( 8 ) 


We can compare the integrands related to various con¬ 
tributions to leTPE as summarized in Fig. 0] The plot is 
prepared using the standard dipole approximation 


Gdip(9') = 


A2 


g2+A2 


(9) 


(with A^ = 0.71 GeV^), which is good to estimate the 
size of various contributions, but is not accurate enough 
for real calculations. 



q(GeV) 


FIG. 4: Integrands for the individual contributions to IbTpe'- 
for the Friar term ([2]) (the black line), the k term (the cyan 
dot-dot-dashed line), the EF term (the red dashed line), the 
Ml term (the green dotted line), and the M2 term (blue dot- 
dashed line). 


III. A PRELIMINARY CONSIDERATION 

As we see from the plot, the Friar term is indeed the 
dominant contribution. However, that is not only the 
value of the contribution that matters. The uncertainty is 


not necessary related to the size of the contribution. The 
Friar term is not only the dominant term for the eTPE 
contribution, but it is also the part of eTPE, which is the 
most sensitive to the parametrization of the form factors. 
The sensitivity of this term to the parametrization was 
discussed in 

The situation with accuracy of the data for the in¬ 
tegration is presented in Fig. [S] taken from [12]. The 
estimations are made within the dipole parametrization. 



FIG. 5: Top: Fractional contributions to the integrand of 
/ft in m as a function of q/A as follows from the dipole 
model. The red dot-dashed line is for the subtraction term 
with G/O). The dashed line is for the subtraction term 
with 1 and the blue solid line is for the term, which 
represents the data. Bottom: the integrand (in arbitrary 
units) is presented with a dashed line and the integral 
Jq dq/q'^ [{Gsiq^))^ — 1 — 2G'e{ 0) q^\ (in fractional units) in 
presented as a function of q with a solid line. The shadow belt 
is for the area in the q space from which 90% of the integral 
Ift come. See [T^j for detail. 

The gray area is for the area, where 90% of the integral 
comes from. One sees that below 0.6A ~ 0.5 GeV there 
is a massive cancellation between the data term (with 
Ge( 9^)) and the subtractions terms (with 1 and with 
G'(0)). The data on Gsiq^) roughly have accuracy at 
the 1% level and below this point the accuracy of the in¬ 
tegrand becomes worse than 1% and increases very fast. 
Meanwhile, above 0.8A ~ 0.65 GeV, the contribution of 
the data themselves (i.e. of the G_e( 9^) term) is negligi¬ 
ble, while the subtraction term with G^(0) dominates. 

Therefore, the integrand at low momentum is deter¬ 
mined by our knowledge of the curvature of (G_e( 9^))^, 
while at high momentum it is determined by the slope of 
(Ge( 9^))^- Both values cannot be measured experimen¬ 
tally, but can be extracted on base of a certain fitting 
procedure. Except for the k term, the integrand of which 
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does not involve the data on the form factors, the other 
data are also sensitive to the fittings. At low momentum 
all the integrands are determined by the slope of the re¬ 
lated terms and thus are proportional to the electric 
and/or magnetic radius squared. Indeed, any fitting pro¬ 
duces a certain systematic uncertainty. 

Unfortunately, the determination of neither electric 
nor magnetic radius is free of questions. Different ex¬ 
tractions from different kinds of data lead to inconsistent 
results. We have already mentioned a strong discrepancy 
between muonic-hydrogen [H and the scattering val¬ 

ues (see, e.g., [a-S|)- The results on Rm are not in a 
perfect shape as well (see, e.g., [Tl|, H^). However, the 
discrepancy in determination of Rm from MAMI Q and 
global data @ is not that important. The most impor¬ 
tant is that the determination of the charge radius Re 
from the scattering and muonic hydrogen is not in agree¬ 
ment, which compromises the whole procedure. Mean¬ 
while, the magnetic radius is even more sensitive to var¬ 
ious systematic effects [iiili. 

As we mentioned, we discuss a rather small correction. 
We anticipate that a different treatment of the higher- 
order proton-size contribution will somewhat shift the 
value oi Re from muonic hydrogen, but not significantly 
(cf. [H). So, we expect that the discrepancy will remain 
and we have to find a way how to use the scattering data 
without a contradiction to the finally extracted value. 

We intend to use fits, rather than the data. Generally 
speaking, to find an appropriate fit is a tricky issue. A 
pure empiric fit is the best suited to fit the data, how¬ 
ever, it ignores various theoretical constraints and con¬ 
straints which follow from some other scattering channels 
(through theory). On the contrary, a fit, which follows 
all those constraint, either may disagree with the data 
or may involve so complicated fit function that it would 
produce instabilities. Therefore, we prefer to use empiric 
fits, but we intend to apply them in a specific way. 

There are two ways to apply empiric fits. E.g., when 
one evaluates the hadronic vacuum polarization contri¬ 
bution to the anomalous magnetic moment of muon (see, 
e.g., [12), there is no subtraction involved. The fit is 
necessary only as a smooth approximating curve which 
is consistent with all the data involved. Its parameters 
by themselves are irrelevant. The other story is when the 
parameters of the fit are really applied. In the integral 
in m they are. The very shape of an empiric fit has 
no physical sense. The parameters must be applied with 
a caution. A scenario, in which the result for the Friar 
correction does not depend on the low-q^ behavior of the 
fit, was developed in [l^ . 

The most plausible scenario to explain the discrepancy 
between the scattering and muonic-hydrogen values of 
Re IS that the scattering data are more or less correct, 
but the accuracy in determination of the radius is overes¬ 
timated. The alternatives, such as a suggestion that the 
data are not correct or the muonic value is wrong, are not 
an option here. If the muonic value is wrong, there is no 
sense in a calculation of small corrections while the dis¬ 


crepancy is that large. The scattering data from different 
experiments are consistent and it is hard to imagine that 
they are completely wrong. 

The integrand involves massive cancellations at low q 
domain, while for high q it is dominated by the contri¬ 
bution, proportional to 

0 )^( 0 ) = ( 10 ) 


If one calculates the integral using scattering data, as 


it was done, e.g., in [13|, the value of G'^(O) is indeed 


taken from there and it is inconsistent with Re from 
the Lamb-shift. If one calculated the whole integral as a 
single thing, then there is no other choice because if we 
take a value of G^(0) inconsistent with the applied fit for 
G_e( 9^), the integral is infrared divergent. 

Following [l2 , our treatment of the Friar term is based 
on two important steps. First, following Hill, we split 
the integration into two parts 


1 = 


dq... = /< 


I> = 


dq... + / dq... . (II) 


Those two terms are to be treated differently. 

The separation disentangles the subtraction at low q 
and the asymptotic behavior at high q, which both in¬ 
volve G^(0), and, in particular, that allows to use for 
the high-g asymptotics the physical value oi Re- This 
value is not fixed, being an adjustable parameter to be 
eventually found from a comparison of the theory and 
experiment on the muonic-hydrogen Lamb shift. In this 
area we integrate the G^ term and the subtractions sep¬ 
arately and while integrating G^ we rely on the empiric 
fits. 

Meanwhile, at the low-q part of the integral we rely on 
the expansion of the form factor 


{GEiq^)f c^l-^q^ + Kq'^ . ( 12 ) 

Even a very bold estimation of the K coefficient allows 
at lo w q a better accuracy than the data by themselves 
(see [l2 and consideration below). 

The final result takes the form 


I = GeiRI + T. (13) 

Such a form of presentation is crucial, allowing a self- 
consistent consideration. As the result, the theoretical 
expression for the Lamb shift in muonic hydrogen keeps 
the desirable form as in o and the higher-order proton- 
finite-size term m contributes to both coefficients ci and 
C 2 . Meantime, we use the empiric fits only in the area 
of the high -5 part of the integrand and the result only 
marginally depends on the details of their extrapolation 
to 5 ^ = 0 (see [12 and consideration below). As a final 
important step, we have a free parameter go which we are 
to vary to minimize the uncertainty of the calculations. 
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IV. CALCULATION OF THE LEADING 
PROTON-FINITE-SIZE TPE CONTRIBUTION 
(THE FRIAR TERM) 


momentum contribution and its uncertainty 

/Fr> = 10 (1 ± 1) ^ . 


(15) 


Let’s consider evaluation of the leading proton-finite- 
size TPE contribution /pr in more details. To deal with 
a lower-momentum contribution 


^Fr< 


dq 

Jo 


{GE{q^)y -1-2 G’e{0) 


we expand the form factor 


{GEiq^)y ^ 1 - ^q^ + ± 1) <z^ (14) 

and, following [H , we estimate the q'^ coefficient with 
help of the standard dipole model (AT'^’p = 10/A'* ~ 
19.8 GeV“^). The dipole formula is not that bad for 
bold estimations. E.g., the G'(0) terms are different from 
the value following from the muonic Re by roughly 7% 
and from those following from hydrogen-and-deuterium 
spectroscopy’s Re and from scattering’s Re hy about 
17%. The uncertainty of 100% for the curvature seems 
reasonable. The value is also consistent with the curva¬ 
ture of various empiric fits (see, e.g., [ssl - l^ : see also a 
discussion below). 

The physical condition on the expansion is qo < 2m-,^ ~ 
0.28 GeV. We are to vary q^ to minimize the final uncer¬ 
tainty and we will check this constraint after we specify 
the optimal value of go- 

The expansion da allows to easily find the low- 


For the low-g part we have exactly followed [l2|, while 
for the high-momentum part we prefer a certain alterna¬ 
tion, which would be more suitable for the calculation of 
the recoil corrections. The higher-momentum part con¬ 
tains the data-related and radius-related contributions 



Gombining the low-g^ and high-g^ contributions, we 
find 


Jpr = 


dq 


lo q 

Jqo r 


{GE{q^)f -I-2 G'e{ 0) q^ 


{GE{q^)f-l 


C'_R|:Pr 


1 

Wo ’ 


(17) 


which is necessary for the decomposition according to 

m- 

To calculate the central value of the data-related term 
in (HD), the integration in [l^ was done over the fits [s^ - 
113 (see Appendix^]) and as the result the related median 
value was taken. Here, we use the most recent empiric 
fit from [13 


GE{q^) 


1 -b 2.909 66Tp - 1.115 422 29 t 2 -b 3.866171 x IQ-^r^ 

1 -b 14.518 7212rp -b 40.883 33r2 -b 99.999 998t^ + 4.579 x lO-^r* -b 10.358 0447r5 ’ 


(18) 


where Tp is defined above in ([S]). We estimate our knowl- comes from a relatively small area around the low limit 
edge of the form factors as at the level of 1% and expect (cf. 0). So, we expect the whole G^ contributions 
that the dominant contribution to the integral for G^ known at the level of 2% and thus we arrive at 


[ $ (<^s(g^))^ X (1±0.02) - ^-b , (19) 

Jqo q OIJq O qo 


where GE{q^) is taken from (fT51) . 


While choosing the fits, we have looked for those which 
produce a smooth function for the all space-like q^ and 
have a good value of the reduced We prefer ‘simple’ 
empiric fits. We believe that it is tricky to determine the 
accuracy for the extrapolation to low q^ and for the sub¬ 
sequent differentiation (to find Re) and, most probably, 
their inconsistency with the muonic value of the proton 
radius is due to an overestimation of their accuracy. Nev¬ 
ertheless, the form factor by itself (not its derivative) is 
fitted by those fits well. The accuracy of the application 


of the fits was estimated in [T3 separately as 

= Wo Ge\^^ (Gdip(9o)) , (20) 

where we expect that 6GE{qo)/GE{qo) — 1%. 

Here, we apply a somewhat different estimation of the 
uncertainty. It is estimated as 2% of the G^ contribution 
into /pr>. The dominant part of the G^ integral comes 
from the low end of the integration area and as far as the 
low-end area dominates the estimations are close. 

Above we explained that it is not consistent to use any 
fits. However, we need to distinguish using of the central 
value of the fits with a pre-estimated uncertainty, as we 
do, and the use of the parameters of the fits to find, e.g.. 
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some derivatives etc. If the fit over a large sample of the 
data has the reduced about unity, that means that 
the fit is consistent with the data at the level of their 
uncertainty, which is about 1%. Such use of the fit does 
not assume any model. Once one intends to take advan¬ 
tages of the large number of data points and to obtain a 
certain value, e.g., G"(0) with a high precision because of 
statistical average, the model comes out and the system¬ 
atic uncertainty due to its application is to be estimated. 
We prefer to stay on the safe side and introduce the un¬ 
certainty at the level of the accuracy of the data which 
grants model independence while applying any fit with a 
good value. 

The scatter due to choice of the parametrization (i.e. 
due to choice of the fits is smaller than the un¬ 

certainty we introduced. The reason for the bigger value 
of the uncertainty is that all the fits basically fitted the 
same data, which were only partly updated from a pub¬ 
lication to a publication. The scatter does not reflect 
the accuracy of the data, but rather the importance of 
their updates. The uncertainties of /Fr< and /Fr> are 
considered as independent. 


Fit 

Ref. 

Type 

Re [fm] 

K [GeV-^j 

(EB 

M 

chain fraction 

0.90 

34.3 


[3^ 

chain fraction 

0.90 

35.3 

(IXSI) 

[33] 

Fade approximation (g^) 

0.86 

28.0 

(IXtI) 

[sa 

Fade approximation (g^) 

0.88 

31.1 

(1X91) 

\m 

Fade approximation (g^) 

0.87 

28.2 

na 

[37] 

Fade approximation (g^) 

0.88 

31.3 


TABLE I: The low-momentum expansion of the fits for the 
electric form factor of the proton applied in [T2|. The values 
are given for central values of the fits without any uncertainty. 
Here: {Gsiq^))^ = 1 - R%q^/3 + Kq^ + .... The related 
values for the standard dipole fit are Re = 0.811 fm and 
K — 19.8 GeV”**. The spread of the central values of the 
charge radius (from 0.86 to 0.90 fm) for different fits (from 
different sets of the data) is comparable with the spread of 
central values of the updated results on the proton radius, 
obtained by different methods, the central values of which 
are 0.84 fm (/rH Qj), 0.88 fm (H&D [J), 0.895 fm (global data 
without MAMI 0— and the results of the other evaluations 
are similar) and 0.88 fm (MAMI i). All the values, but from 
muonic hydrogen, have uncertainty at the level of 1%. The 
muonic one is much more accurate. 


as we mentioned, it simplifies a comparison with [l3l| . 
Note, that is not an improvement, but a simplification. 
The reason to check the median using various fits (some 
of which are partly out of date) was to check that the 
result does not much depend on low-g^ behavior. The 
spread of the values of the radius was sufficient to prove 
such model independence, which was also confirmed by 
the use a fit from [s^, which, despite of unreasonable 
low-q behavior, produces a quite reasonable result. Only 
after such tests are performed, the use of the most recent 
fit is possible. 

The extracted value of the proton radius depends on 
the sep aration parameter qq in almost the same way as 
for [T^. The procedure (see [l^ and sections below) al¬ 
lows to vary the value of the separation parameter qq and 
to find an optimal one. The results with different values 
of go are consistent. The deviation from the former treat¬ 
ment of the eTPE [l^ is large enough (at the level of the 
uncertainty) to be considered seriously. The dependence 
Rsiqo) is shown in Fig. El 


flE(fm) 



FIG. 6: The value of the proton charge radius Re extracted 
[T^ from measurement [J] of the muonic-hydrogen Lamb shift. 
The plot follows [l^ . The horizontal belt is for the original 
value of the parameter. The other figure [of a more compli¬ 
cated shape] is for the radius determined in [I^ , presented as 
a function of the separation parameter go. Every area corre¬ 
sponds to one sigma. The internal and darker part of each 
figure is related to the partial uncertainty due to the calcula¬ 
tion Jft, while the whole area is for the total uncertainty of 
the radius. The final result of Re = 0.840 22(56) fm was 
obtained there at the optimal value go — 0.152 GeV. 

Minimizing the uncertainty we arrive at 


To facilitate further numerical evaluations with includ¬ 
ing the recoil effects, we choose here instead of the median 
the fit [s^, which is the most recent one. This modifi¬ 
cation is not important because our estimation of the 
uncertainty is substantially larger than the scatter. The 
most of the integral (for the term) comes from the low 
limit. Here, we choose to assign the 2% uncertainty to the 
whole term. We have checked that this modification 
is of marginal importance for the Friar term, however. 


/p, = - 23.5(3.3) GeV^ , (21) 

3 go 

where the numerical value corresponds to go — 
0.147 GeV. It is important that the new optimal value is 
quite below 2m,r which is the requirement for the validity 
of the expansion (HI. 

While the decomposition term and into the con- 






















stant / is somewhat different^, the related numerical re¬ 
sults at Re = 0.84 fm are in a perfect agreement, being 
17.6(3.4) GeV^ [13 and 17.7(3.3) GeV^ (this work). 
The approach applied here is a modification of the one 
from [l3 and the results for the radius are expected to 
nearly coincide at Re = 0.84 fm. The latter effectively 
means that the extracted radius is to be approximately 
the same (with a marginal deviation). (We widely use in 
this paper the ‘reference’ value of various contributions 
at Re = 0.84 fm—that is the approximate value consis¬ 
tent with extractions from the muonic hydrogen [l|, [T3 
at the level of 10“^ which is completely sufficient taking 
into account the uncertainty of our calculations of the 
Friar and eTPE integrals at the level of 20%.) 

Those numerical values should be compared with 
the results which can be obtained if one directly ap¬ 
plies empiric fits to 0. The related results are 
/ft = 23.1(1.1) GeV“^ (after ap plyi ng empiric fi ts (doll l. 
22.9(1.2) GeV“^ (after applying [IJ] empiric fits [351 and 
[13 — that is similarly to consideration in [13 , which was 
done for the complete Istpe) and /ft = 24.3(7) GeV“^ 
(after applying empiric fits [lOj)- The uncertainty of our 
result is essentially larger than theirs and the the shift 
from the quoted values to ours is larger than the un¬ 
certainty. We believe that the uncertainty of the cited 
results was underestimated and suffered from an incon¬ 
sistency (once we accept the muonic value for the proton 
radius). 

We have to recall the geometrical meaning of the Friar 
term 

-^Fr = ^ y d^r J dV'|r - r'|3p(r)p(r') , (22) 

where in a non-relativistic approximation p{r), being a 
Fourier transform of GF;(q), is the charge density in¬ 
side the proton. The muonic-hydrogen data assume a 
smaller proton, thus the Friar term adjusted to those 
data should also be smaller (see also discussion on that 
in [i3)- The radius is smaller by about 6% than the re¬ 
sults from the empiric fits and the Friar term (the central 
value) is smaller by about 15%, which seems consistent. 

Indeed, in a rigorous consideration, the Fourier trans¬ 
form of Ge is not the charge density in common sense. 
However, the presence of relativistic effects does not af¬ 
fect the fact that the value of the Friar term is determined 
to a certain extent by the proton size, it just makes this 
relation more complicated and involves corrections. The 
non-relativistic part with a clear geometric sense remains 
and it is dominant. 

One may wonder whether splitting the area of the inte¬ 
gration maintains the continuity of the form factor. We 


^ The result of [T3l is l-pi = ~ 22.2(3.4) GeV“®, where 

the numerical value uses the optimal value qo ~ 0.152 GeV. 
(Indeed, since the uncertainty of /> for the Friar term is defined 
here and in [T3l differently, the optimal values, which minimize 
the uncertainty, are also somewhat different.) 


cannot check it a priori because at low q'^ we use the 
actual value of the proton radius. However, once the 
complete calculation is performed, we find the value of 
the radius and can to check the continuity. Such an a 
posteriori test is plotted in Fig. jT] The red area is for 
the low-q^ expansion and the green one is the fit from [s^ 
suggesting the uncertainty of 1%. They are consistent. 


G/Gdip 



0.0 0.1 0.2 0.3 0.4 0.5 


FIG. 7: The continuity test for the proton electric form factor 
normalized against the standard dipole fit. The red area on 
the left with fast increasing (to the right) width is for the low 
momentum approximation, while the belt with approximately 
homogenous width on the right is for the empiric fit [^ we 
used, with the uncertainty of 1%. 

The method suggested in [l3 pretends that the re¬ 
sult does not depend on the parametrization of the fit. 
To check that, a number of fits have been applied (see 
Fig. E]). Those include the chain-fraction fits (in q'^) [^ 
(the blue dashed line), the Pade-approximatioii (in g^) 
fits [13113113 (the green dot-dashed lines) and [I 4 I (the 
green bold solid line), and a fit with the inverse polyno¬ 
mials in q [13 (the solid red line). The fits except for 
[13 have reasonable low q behavior. They are analytic 
in q^ around q^ — 0 and the related value of the charge 
radius varies from 0.86 to 0.90 fm, so the spread is com¬ 
parable with the distance to the muonic value of 0.84 fm 
(see Table [J). The high-g^ asymptotics is proportional 
to g“^ for the chain fractions and to g“'* for the Fade 
approximations. Therefore, the parameters and the gen¬ 
eral behavior (far from the characteristic go area) vary 
very much. Nevertheless, the spread of the values for the 
integra l J^ '> is smaller than the estimation of the uncer¬ 
tainty [ 13 . The scatter of the results with inclusion of 
the fit with unrealistic behavior [13 is somewhat larger 
than the uncertainty but still acceptable [13 • 

V. CALCULATION OF THE RECOIL 
CORRECTIONS 

After discussing the leading eTPE term (within the ex¬ 
ternal field approximation), we are prepared to consider 
the eTPE recoil corrections to it. We start from ([7]). As 
we see from Fig. 21 the dominant proton-finite-size terms 
are the EF and the k ones. The former is the proton- 
size term, while the latter is not related to the proton size 
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FIG. 8: The fractional deviation of the fits (see in App. 
for details) and the data for the electric form factor Gsiq^) 
from the standard dipole form factor. The data are taken 
from [ 3 ^ , where the proton form factors were extracted from 
experimental data on the elastic electron-proton scattering. 
The fits include the chain-fraction fits (HU and (IA3I) (in q^) 
[ 35 I (the blue dashed line), the Pad e-app roximation (in q^) 
fits dMll, ISl), iTMll . and (Hi) p, IM a (th e green dot- 
dashed lines) with the most recent one from shown with 
a green bold solid line. An exceptional fit (IA12II with an 
inverse polynomials in q [ 3 ^ is shown with a solid red line. 


directly. Domination of the EF term is due to q/M ex¬ 
pansion. Comparing the coefficient functions at the area 
of m <C q M we find 


== (23) 


which explains why the EF term is dominating. 


An important area of the integration is related to g ^ 
m and it is useful to present a more detailed asymptotics 
at TO <C M, q M (exactly in m/q) 


/ef 


^\/4to^ -I- ( 7 ^ — + 8 to"^ + q^ — 2m^ — 2rn^q^ 


(24) 


Decomposing the EF term according m we find 


and in a close analytic form 


Ief= 


® dq 

4 


+ 


0 9 

90 9^ 

Qo 


{GEiq^)Y -l-2G'EiO)q^ 


[GE{q^))^ - 1 /ef(9^) , 


ChI:EF--11 ^/EF(g"). 


(25) 


That is very similar to the expression (HZJ for the Friar 
term. The evaluation of the integrals is done in a similar 
way and we apply the same effective description for Ge- 
The only difference is that in former case all the integrals 
except for the G^ term were calculated analytically, while 
here we have to do that numerically. 

The other important term is the k term. It is a 
proton-structure correction, but not a proton-finite-size 
one. Since we decompose elastic TPE into the point¬ 
like Salpeter term and a proton-structure term, the sub¬ 
traction in the latter is due to a proton without any 
anomalous magnetic moment. Therefore, there is a con¬ 
tribution to eTPE which does not vanish in the limit 
Re = Rm = 0, which is the k term. It does not involve 
the parametrization of the form factors and may be found 
with a high accuracy. We obtained the result from both 
numerically 


/« =-2.850 23 GeV^ (26) 


h. = 


M 

-3(1+ k) In — 
m 


16 

+ (!-«:) +B- (^) 


(27) 


where the value of B ^ 0{1). Our result for this value is 


B 


k-2, 2M 2 +11k 

-In-^- 

4 TO 48 

'k-3^^2M 43(1-bAt) 

16 TO 384 



As one sees (cf. (I^5)l ). this term is ~ M~^ in contrast 
to the leading part of the EF term, which is ^ M~^. 
However, the non-vanishing value of k, which enters the 
combinations such as k{1 + k), produces a certain en¬ 
hancement. 

Let’s now consider the magnetic finite-size terms, /mi 
and Im 2 - At low-g we use expansion 

(Gm( 9"))" = Mp (^1 - ^9'[1 ± 0.4] + K<^'Pq\l ± 1)^ , 

(28) 

where /ip = 1 -|- Kp ~ 2.792 85 is the proton magnetic 
moment in the units of the nuclear magnetons, for the 
magnetic form factor, similar to that for the electric one. 
However, there is a difference. While for the electric form 
factor the value of Re is the adjustable parameter for the 
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extraction of the charge radius from muonic hydrogen, 
the rms magnetic radius Rm should be a constraint. We 
use here the dipole value of the radius 

, 12 
jdZ _ _ 

— ^2 ’ 


but with a conservative uncertainty of 40% (to R\j)- 
As for the high-g, we apply the fits from [s^, which 
gives not only the fit for the electric form factor m, 
but also a fit for the magnetic one. The latter is of the 
form 




1 - 1.435 73Tp + 1.190 520 66t^ + 0.254 558 41t 3 
1 + 9.707 036 8lTp + 3.7357 x 10-V2 + 6.0 x IQ-SrS + 9.952 7277t* + 12.797 7739r5 ’ 


(29) 


go [GeV] 

[GeV-i] 

/eTPE [GeV-3] 

/eTPE [GeV-3] 

0.1 

3.362 20 

-44.52 

16.4(6.0) 

0.2 

1.714 65 

-12.79 

18.3(4.0) 

0.3 

1.175 25 

-2.597 

18.7(5.8) 

0.4 

0.91153 

2.527 

19.1(7.6) 

0.5 

0.757 09 

5.790 

19.5(9.2) 


TABLE II: The results I and IbTpe for fits (I18II and (I29II at 
various go- The last column is for the numerical value of Rtpe 
at Re = 0.84 fm. 

The decomposition into I and and further numer¬ 
ical evaluations of those values are similar to those for the 
EF term. 

The numerical results for eTPE as a function of differ¬ 
ent values of the separation parameter are summarized 
in Table HIl The minimization of the uncertainty leads to 
go — 0.147 GeV as the optimal value. 

The result is 

4tpe = [-24.14(3.27) -h 59.31 r^] GeV^ , (30) 

where the error of 2.90 GeV“^ corresponds to /eTPE<, 
while the error of 1.53 GeV“^ relates to /eTPE>- 

The individual contributions to the eTPE correction 
are presented in Table Hill while the individual contribu¬ 
tions to the uncertainty are summarized in Fig. |9l We 
see that the individual electric and magnetic recoil contri¬ 
butions are at the level of 10-15%, however, they have a 
trend for a strong cancellation. The cancellation is rather 
accidental and depends on the value of the anomalous 
magnetic moment of the proton Kp, the mass ratio m/M 
and the proton radius. The terms differently depend on 
the parameters and the cancellation of the central values 
does not produce a massive cancellation of the uncertain¬ 
ties. 

However, the latter are much smaller, than the level 
of the 10% of the Friar-term uncertainty. The size of 
the contribution is not that important as its sensitivity 
to the parametrization. It is also important, that any 
uncertainties at the level of 10% of the Friar-term uncer¬ 
tainty would be important only if it is correlated to the 
main uncertainty and we have to sum them up linearly 


rather than as the rms sum. That leaves us with the EF 
term as the main source of such an uncertainty since it 


Gontribution 

Cnl [GeV-i] 

/eTPE [GeV 

/eTPE [GeV 

Fr 

2.271 

-23.5(3.3) 

17.68(3.31) 

K 

0 

-2.850 

-2.850 

EF 

0.039 

1.26(4) 

1.97(4) 

Ml 

0 

1.76(3) 

1.76(3) 

M2 

-0.0006 

-0.838(4) 

-0.849(4) 

eTPE 

2.309 

-24.1(3.3) 

17.71(3.27) 


TABLE III: Individual contributions to the Etpe at optimal 
go — 0.147 GeV. The last column is for numerical value of 
7eTPE at Re = 0.84 fm. 


is expressed in terms of the same parameters as the lead¬ 
ing term. The contribution to the integral by itself can 
come from a various area of g, however, the area sensitive 
to the parametrization is around go. The expansion at 
low g (< go) is the most uncertain at the largest ‘low g’ 
available, while the uncertainty of the high-Q area (> go) 
is mostly determined by the smallest ‘high g’ available. 
The contribution comes from a more broad area than that 
around g ~ go ~ 0.147 GeV. In the limit m <C go “C M, 
the suppression of EF behaves as qo/M with small coeffi¬ 
cients (< 1) as follows from the asymptotic behavior (1231) . 
The optimal go is found to be go — 0.147 GeV and thus 
go < That leads to a suppression of the EF term 

in respect to the Fr term by a factor of m/M (^ 10%) 
with small coefficients (< 1) partly because of additional 
softening of the integrands by factors g/2m^ < 1 and 
softening of the character of the g-integration. 

As an important test we check that the scatter due 
to use of the empiric fits (except for those from Hi) is 
below the uncertainty. The fits from [li shift the result 
somewhat farther, but the value is acceptable because of 
the unreasonable behavior of the fits at low g. 


VI. CONCLUDING REMARKS 

The extraction equation is of the similar form as for 
the Friar term in [1^ 

















0.2 


0.3 


0.4 


q(GeV) 


0.5 


FIG. 9: Individual contributions to the uncertainty J/stpe 
in relative units of the dipole value for various go- The total 
error SI is presented by the thick solid line, other types of lines 
are described in the caption of Fig. [4] The recoil contributions 
are very small and to put all the contributions into the same 
picture we use different normalization. The values of 57eTPE 
and SIft are multiplied by 0.1, which allows to use a linear 
scale (instead of the logarithmic one). The recoil corrections 
are [anti] correlated with the Friar term, which somewhat 
reduces the total uncertainty. 


Fit 

Ref. 

Type 

1 

> 

7eTPE 

[GeV^^] 

(TATli. (TA^ 

m 

chain fraction 

-24.86 

16.99 

(IMt. (TAl 

[3^ 

chain fraction 

-25.01 

16.83 

(lA^. (jAeJ 

[33] 

Fade approx, (q^) 

-23.66 

18.18 

(lATli. (lA^ 

m 

Fade approx, (q^) 

-24.18 

17.66 

((A9li. (IaToIi 

[3^ 

Fade approx, (q^) 

-23.86 

17.98 

(TMli. (IMTl) 

m 

Fade approx, (q^) 

-23.87 

17.98 

(0, (I29j 

[37] 

Fade approx, (q^) 

-24.14 

17.71 

(iMli. (IAT 3 I) 

m 

inv. polynom. (g) 

-29.31 

12.54 


TABLE IV: Numerical results for central values of I and /stpe 
at the “optimal” qo — 0.147 GeV from various fits. The last 
column is for numerical value of Ibtpe at Re = 0.84 fm. 




( 0),2 _ ‘2.4:{Za)mr 


1 


1 - 


2A{ZoL)mT 


c 


r%{Qo) 


C'_R|(9o)(^bV + (?eTPE(go) 


~ ^eTPE(90 


(31) 


where /g-ppE i® original result for leTPE and leTPE 
is the one from which we removed the term. The 
coefficient Cr 2 is due to a small ‘renormalization’ of the 
RY; term. 

In the case of the Friar term the ‘renormalization’ is 

, , 8(Za) rur 

1 -)> 1-^- - — 

TT go 

and the recoil effects weakly affect it. For /gppE we use 
here the result from [I^ with the Fl term restored as 


suggested in (l3 | 

= (^27.5(1.6) - 4 . 5 ) GeV“^ 

= 23.0(1.6) GeV"^ . (32) 

That is related to the evaluation in [l[. 

Eventually, we arrive at 

i?£; = 0.840 22(53) fm , (33) 

which is to be compared with the original result Re = 
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0.840 87(39) fm from [I| (which we use here as the ref¬ 
erence value i?®) and the result Re = 0.840 22(56) fm 
corrected for the self-consistent evaluation of the Friar 
term [H- Note that in both papers on the self-consistent 
reevaluation (0 and this paper) only a correction was 
found, while all the other theory remains the same as 
used in [I|. 

Various extractions of the proton charge radius are 
summarized in Table El The main effect of the self- 
consistent evaluation comes from the correction for the 
Friar term. The results of [l^ and this paper, while cor¬ 
recting the Friar term, have a marginal difference which 
confirms that the assumption, that while calculating the 
term of J> a small area of the data around the low 
limit dominates, is valid. The result of the complete eval¬ 
uation is marginally different from the result of the re- 
evaluation of the Friar term only. The strong cancellation 
between corrections from various contributions is rather 
accidental, while the marginal shift of the uncertainty 
is rather expected, because the corrections to /etpe are 
small and less sensitive. 


Re 

/eTPE 

Comments 

[fm] 

[GeV-3] 


0.840 87(39) 

23.0(1.6) 

Re from [^, /exPE from 

0.840 22(56) 

17.7(3.4) 

Ipr from [^, free from 

0.840 21(53) 

17.7(3.3) 

Ipr from this work, Aec from 

0.840 22(53) 

17.7(3.3) 

/eTPE from this work 


TABLE V: The results of various evaluations of the proton 
charge radius. The theory follows except of the treatment 
of the eTPE term. The eTPE term in [![ was found following 
M- The result from [a dealt with the correction of the 
Friar term only. The recoil part of the eTPE is added to the 
Fr contributions accordingly to [a.0. One of the result 
of this paper was obtained also by correcting the Fr term 
only, while the other is a result of the complete results of the 
self-consistent evaluation here. 


The uncertainty of the result (1331) comes from our 
knowledge of the curvature of Gb)^^) at low q and on our 
accuracy of the data on Gb around g = 0.1 — 0.2 GeV. 
Any improvement of those will immediately improve ac¬ 
curacy of determination of the proton charge radius from 
muonic hydrogen. 


in part by DFG (under grant # HA 1457/9-1) and by 
RFBR (under grant # 14-02-91332). 


Appendix A: Fits for the electric and magnetic form 
factors of the proton applied in the paper 


The results in the paper have been obtained by evalu¬ 
ation with the fits from [3^ . However, it is important to 
be sure that the result does not depend on details of the 
fits and in particular on details of their expansion around 
q = 0 and their asymptotics at high q. Various fits for 
Ge and Gm are applied in this paper for test purposes. 

The fits include two pairs of the chain-fraction fits. 
Both are from Arrington and Sick, 2007, [s^. The first 
one is completely based on [sj 


G_e(9^) 


1 


3.440=^ 


Q.178Q^ 

1 .212Q^ 

I I 1-176Q^ 
1-0.284Q^ 


(Al) 


and 


Gm(Q^) 

Fp 


1 


1-b 


3.173Q^ 



0.314Q^ 

1- 

1.165Q^ 

- 7 ^ - 


l-1.087Q^ 


(A2) 


where, Q is the numerical value for the momentum trans¬ 
fer q in GeV. 

The other pair is obtained in 
photon correction according to 


M I by applying the two- 

41| 


and 


Gsiq'^) = 


3.4780^ 

Q.140Q^ 

1.128Q^ 

1-0.233Q^ 


(A3) 


GM(g^) 


3.224Q2 

0.313Q^ 

^ 0.868Q^ 

, , 4.278Q^ 

1-1.102Q^ 


(A4) 


Other fits are with the Fade approximation in orig¬ 
inally introduced by Kelly, 2004, such as [s^ 
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^ 1 + 10.98rp-b 12.82r2-b 0.863r3 ^ ’ 

and 

GM(.q^) _ _1 -b 0.12rp _ 

tip “ 1 -b 10.97rp + 18.86r2 -b 6.55r3 ' ^ ^ 

The Pade approximations have been later developed 
by Arrington et ah, 2007, [s^ 
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Ge- 


1 H” 3.439Tn — 


1.602t2 + 0.068t3 


1 + 15.055Tp + 48.061r2 


99.304t 3 + 0.012r4 + 8.650 t5 


(A7) 


and 

GM(g^) ^ 

In paper 
presented 

GE{q^) = 


1 - 1.465tp + 1.260r2 + 0.262 t3 


1 + 9.627r„ + 11.179r4 + 13.245 t5 


(A8) 


Ip I I 

by Alberico et al., 2009, one electric fit is 


1 - 0.19t, 


p 


1 + 11.12rp + 15.16Tp2 + 21.25Tp3 


and two magnetic ones 
Gniq^) 


1 + 1.53t„ 


Mp 

GMjq^) 

Mp 


1 + 12.87tp + 29.16r2 + 41.40r3 


(A9) 


(AlO) 


1 + 1.09r„ 


1 + 12.3lTp + 25.57rp2 + 30.61 t3 


GM(g^) 

Mp 


0 ) and 
3711 and 


The most recent Pade-approximation fits (see 
are presented by and Venkat et al., 2011, 
have been used in the main body of the paper for all the 
evaluations. 


One more pair of the fits with is with the inverse poly¬ 
nomials in q from Posted, 1995, 


GE{q^) = 


1 


(All) and 


1 -f 0.62g -t 0.68Q2 + 2.8Q3 + 0.83Q4 

(A12) 


= (l -f 0.35Q -f 2.44Q2 + o.SOQ^ -f l.OdQ^ -k 0.34Q3) 


-1 


(A13) 


That is a phenomenological fit designed to be used for 
medium and high q. It is not expected to be appropriate 
at low q. Providing a reasonably good approximation 


at medium momentum transfer, the fit apparently has 
incorrect low-g behavior and incorrect analytic properties 
such as a branch point at q^ = 0. 
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